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SUMEARY 


In order to determine the effect of a non-linear 
spring characteristic on the vibrations of a torsional 
system, investigations were made to determine the response 
curves of a given system when the drive was made through 
am non-linear coupling and when made through a comparable 
shaft. The results showed the maximum amplitude for the 
former cas6é to be about one-fourth of that for the latter. 
However, when the difference in damping was acounted for, 
the non-linear coupling showed a reduction in amplitude 
of vibration to about one-half that of the comparable 
linear system. An analytical solution inelucing danping 


is given, and suggestions are made for further study. 
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While couplings having non-linear torque-deflec- 
tion characteristics have been used for some time in in- 
stallations where it is desired to reduce torsional vi- 
brations, and while the accomplishment of this mission 
is generally conceded, there is no data avallable for 
evaluating their effectiveness in thie regard. It is de- 
sirable, when considering an installation of this type, 
to know not only if an improvement will be realized, but 
also how much improvement. It was, therefore, with this 
point in mind that the subject investigation was undeér- 


taken. 























— ammeter (load amps.) 
~— pedestal bearing 
Yhitworth drive 


S 
é 
SC SF GS > 


— driving dise (referred to coupling), graduated 
0-360 degrees 


Cc - @€ = non-linear coupling 


d — driven dise (referred to coupling), graduated 
0-560 degrees 


GC generator (GeCe, 115 ve, 1/20 HePe at 1725 ReP We) 
IL, —- lampbank load 


Mi — motor (GeCe, 115 Ves 2045 AMPSBe, 1/4 HeP. at 
1725 ReP ole) 


P — breaker points, attached to graduated dial 
R — rheostats for speed controeal 

r — rheostat for control of generator field 
Ra = mounting rail 

S — stroboscope 

T -—« torsiograoh | 

Vew- voltmeter (generator field voltace) 

Vew= voltmeter (generator terminal voltage) 


X w= “stroebctac” {not shown) 


Power: G Ve GoeCofor satroboscope and torsiograph 
115 v. d.c.—for motor drive and generator field 


115 ve. €0 cycle aeCemfor strobotac 








The Whitworth drive is an inversion of the slider 
crank mechanism designed to impress an oscillatory torque 
on dise "b". The offset between the shaft to "a" and the 
shaft from "b" was kept constant at 1/8" throughout. 

The coupling used is shown in detail in Figures 3 
and 4. It may he seen that the segments holding the spring 
strips are cut with first a plane portion and then a ecylin- 
Grical portion. The radius of curvature of this latter part 
is 6 inches. The non-linearity is produced by the wrapping 
of the spring stripa upon the cylindrical segments thus: 


This shortens the unsupported length of the strip and s0 
stiffens the coupling. The spring strips have a free slid- 
ing fit in the straight portion of the segments in order to 
allow for the axial motion necessary when the coupling is 
deflected. Clock-spring steel .025 inches thick and 3/8 
finches wide was used for the spring strips. The “coupling 
lengths" referred to in the discussion to follew are the 
unsupported lengths of the spring strips when not deflect- 
ed. This length was varied for the different runs by slid- 
ing the assembly "G", "ad", and one coupling flange "c" to- 


ward the motor end of the set-up. The ends of the strips 
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were then cut off to allow only the slight projection beyond 
the coupling flanges shown in Pigure 2. Parts "™", “B", and 
"os" were bolted rigidly to the mounting rail "Ra". For the 
runs made with ordinary linear shafts only the soring 

strips were removed from the set-up as shown in Figure l. 
The coupling flanges “c" were allowed to remain in place 
and the shaft connected between them by means of chucks. 

Two intermediate bearings were used to support the shafts 
tested. 

The shaft of motor “" projects through the out- 
board bearing to within a short distance of dise "P". A eam 
on this shaft operates a set of breaker points mounted on 
"p", breaking the cireuit through these points once each 
revolution. ‘The primary of the stroboscope "S" is in series 
with this circuit, causing the stroboscope to flash each 
time the circuit is broken. By rotation of disc "P" to the 
desired position this flash may de made to occur at any 
position of the shaft during a revolution. A knurled lock- 
ing nut is used to lock dise "P" in the position desired. 
The effect of thie mechanism is to sallow simultaneous read- 
ings to be taken of the position of dises “b" and “d" at 
various points throughout a revolution of "b", these points 
being determined by the setting of dial "Pp" which is gradu- 
ated. ‘Thus a scan may be made of these instantaneous posi- 


tions throughout a revolution. A cirenuit from the breaker 
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points also operates the dead center timing pen of the 
torsiograph.e 

The torsiograph is driven by ea belt from the shaft- 
ing between "b" and "B". In order to reduce slippage this 
section of shafting was wound with friction tape. In this 
condition the shaft diameter was 11/16 inches. ‘hen run- 
ning the linear shafts it was found that the amplitudes of 
vibration were too large to be recorded by the torsiograoh 
as installed, even at its lowest amplification. This was 
remedied by driving the torsiograph from the test shaft it~ 
self, which provided the necessary smaller driving diameter. 
This shaft, taped, measured 7/32 inches. 

The generator field was separately excited, the ex~- 
citation being maintained at 100 volts. The lampbank load 
was not changed during the investigation; hence the load in-~ 
creased with the speed of operation. The lamps of the lamp- 
bank were painted black in order that they would emit no 
Light. 

The strobotac was mounted in the position indicated 
by “X" in Figure 1, so as to throw its light on the spring 
strips. “hen ranning the shafts, the strobotac was brought 
to the near aide (in Figure 1) of the mounting rail and its 
light directed on that face of disc "ad" which faces the 
coupling flanges. The atrobotac was used for speed measure] 
ment only, and should not be confused with the stroboscope 


wrich was used in obtaining dial readings as described above. 
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For all runs involving the use of the stroboscope 
or the strobotac, the laboratory was made as dark as pos- 
sible in order to facilitate the taking of readings. 








PROCEDURE 


With the apparatus set up as shown in Figure l, a 
preliminary run was made to determine the speeds at which 
{t would be desirable to take readings. This was done by 
increasing the motor speed by small increments, each time 
"stopping" the rotational motion of the coupling with the 
strobotac. In this state the amplitude of the vibratory 
motion in the coupling could be observed by eye. Soeeds 
were thus chosen which gave desired increments of ampli- 
tude change. The critical speeds were also determined 
during this run. 

The strobotac was then set at the lowest speed at 
which an amplitude determination was desired. The motor 
apeed was adjusted, by means of the motor rheostats, until 
rotational motion of the coupling was "stopped" by the stro- 
botac light. Since the speed of the motor had a tendency to 
vary, it was watched constantly and readings were taken only 
when the speed was steady and at the desired value. The ro- 
tation of the motor shaft caused the stroboscope to flash 
once each revolution, as previously explained. By the 
light of these flashes the dials on the driving and driven 
discs (designated as "b" and “a" in Figure 1) were read 
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and recorded. ‘The eorresponding setting of dial "P"™ was 
also recorded. Dial "P" was then rotated in increments of 
15 degrees, thus causing the stroboscope to flash at suc- 
cessive positions during the cycle. Readings of dials “b” 
and “d" were taken for each position. When a complete cycle 
of readings had been obtained the voltmeter and ammeter read- 
ings were recorded and a record of the motion of disc “b" 
taken on the torsiograph. The speed was then increased to 
the next desired value and the procedure repeated. Headings 
were thus taken at speeds ranging as close to the critical 
spesd as possible. Readings at the critical speed were not 
possible due to the difficulty of maintaining constant speed. 
After passing the critical, the speed was increased to the 
maximum value at which readings were desired, and then ree 
duced in steps, taking readings at each step as before. 
This was continued until the point of instability at de- 
creasing speed was reached. 

The motor was then stopped and a reading taken of 
the relative position of dials "b" and "a" with no twist 
in the coupling. This served as a reference reading in 
determining the actual position of one dise, relative to 
the other, from the readings previously taken. 

Following this a bar was attached to the shaft be- 
tween "BR" and “e" and the end of the generator shaft was 


rigidly clamped to prevent turning, by means of a pedestal 
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secured to the mounting rail. Using a spring scale known 
increments of torque were applied to the shaft and the ree 
sulting angular deflection read from dial "b". This was 
done for increasing and decreasing values of torque to ob- 
tain data for the hysteresis loop of the unit. This com- 
pleted the readings taken for a given coupling length. 

The entire procedure was repeated for each of the 
coupling lengths. 

The procedure for each of the shafta (1.e., “linear 
couplings") tested was identical, with the exception that 
the hysteresis loop was taken before the run in order te 
determine that the shaft stiffness was of the desired mag- 
nitude. This stiffness was also checked after the run to 
assure that no change had occurred due to the stresa of the 
large deflections encountered. 

Due to the small output of the generator and the 
aifficalty of determining its efficiency at this light load, 
a series of runs was made with the motor and generator 
coupled through a short, stiff shaft. The input to the 
motor was determined with no load and no field on the 
generator, and again determined with the load and field 
adjusted to that used in making the coupling runs. ‘The 
difference was assumed to be the load input to the sener-~ 


ator. This difference was cetermined for various speeds 





and the results plotted (see Figure 18). While the accuracy 
of this procedure is not good, it is considered satisfactory 
inasmuch as it contributes only a small fraction of exciting 


torque applied to the driving disc. 
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THEORY 


1. Notation 


component of steady state amplitude, see defini-e 
tion of s and ye 

viscous damping, factor, pound inch seconds per 
radian. 

component of steady state amplitude, see defini- 
tion of s and y. 

the eccentricity between dises “a” and "bd". 
moment of inertia of mass or combination of 
masses indicated by subseript, see Figure 5. 


mass moment of inertia of equivalent system 
Ty1p 
+ 
qt tT, 


torsionsl stiffness of shaft, pound inches per 
radian. 

initial slope of non-linear torque-deflection curve. 
coefficient of the cubed term in the expression 


for the non-linear torque-deflection curve. 


Ke 
3/4 ey 


-S 
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om 


rotational speed in revolutions per minute = Ws 
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i,* natural frequency expressed in revolutions per 


= oO 
minute + Ww, 
“a & 


r * the crank radius of the Yhitworth linkage. 


R * frequency ratio = = r— 
n n 


;*= ye = a” + of 

Tt“, * tT, 

Tp™ the coefficient of the sinusoidal component of 
driving torque applied to dise “b". 

Tp'™ total driving torque applied to dise "yp". 

Tp” damping torque. 


IyIg 
Tg” equivalent torque * ———— 
i, *f 
1 2 
as a load torque. 


T* inertia torcue. 


y * angle of twist in radians. 


4 
jo = “amplification factor" * oe 
J(i-R")® + YR 





= &D 


i 
®,, wn 





y = dimensionsless damping coefficient = 
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2. Derivation of amplitude equation 





In order to reduce the complexity of the analysis 
the system must first be simplified. The steps of this 


simplification may be indicated as follows: 





Figure 5 


The letters indicating the original masses of the system 
above are the same as those employed to designate these 


parts in Figures 1 and éd. 


Sinee the masses "mM" and "a" are rigidly connected, 


= + 
I, Tay I. 


likewise, 


_ I, ie” 


5 c 
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and Ip* I, + Ig * Ig 


Equating the kinetic energies of system 1 and its component 
parts: 
Yet, +*1/2 yeh * V2 Ty Way, 
AV. AV. 


Bat @ = = @ 
Mov. a — 


s+ % 7H, 2% 


we now have a two mass system with moments of inertia I, 
and Io» connected by a spring of stiffness k. 

Let a periodic torque of magnitude T, be applied 
to mass 1. Then it may be shown (Reference 1, page 35) that 
the equation of motion of the two mass system is; 

Iylg + py + uy = Tet, ~~ Ff (wt) 


> 


Ty Tp ty " I, 
Where y * the relative motion between dises “b" and "a" 
(viz., the angle of twist). 
This equation is identical with the equation of absolute 
motion of the system E if 


= Iy1 
:. ite 


I, * tp 


and T," Ipty 
Ty Tp 
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Therefore we may write as the equation of motion of the 
equivalent system: 

Tpy * by + ky = Ty f(t) (1) 
tT, f( wt) may be written as a Fourier series, the motion 
due to each component evaluated, and the net resultant mo- 
tion obtained by superposition. from this point on we there- 
K sin wt, 
which satisfactorily describes this torque for the case 


fore will consider only the fundamental component T 


mnder consideration. 

In the ease of a linear system (k * constant), we 
heave a linear differential equation, which may be solved 
for y, giving the well known expression for maximum ampli- 
tude; 

Ty /k 


Tmax c ‘iar? )? bs 52,2 - (3 (tf) 


However, for the non-linear system |i * r(y)] no 
such simple general solution exists. The following deriva- 
tion is based on an article by Appleton (Reference 2). 

First, it is assumed that the terque deflection 
curve of the torsional spring can be expressed by an equa- 
tion of the form: 

'* ki + ny” 

The ehoice of this form of expression for the spring charec- 


teristic is arbitrary, and is chosen primarily for conveni- 
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fence. However, results of this investigation show it to 
give satisfactory agreement with the actual curves obtained. 

Substitating this expression for our spring torque 
into equation (1) we have: 


Iny * by + kyy * key” = Ty sin wrt, (2) 


£ 
— = .2 
or I,y + by + ky (1 + KY ) = T, sin @t .« 
Dividing by Ips we obtain: 


wy Fi 22 
y* Ty* 7 (1 EY )= TE sin Wet 
Tk 
k 
° z 2 
+ any * my (1 * 


or y’) a 0 yee sin @pt, (5) 


taf 3 


—_ J” Ky : 


Neglecting the transient term, let the trial solution be: 


y = a sin @,t + ¢ cos apt (4) 


r 


then ¥ = @ @pCos Wt - ¢ Wp sin Wet (5) 





and J = - & We sin Wt - © We cos wet (6) 


In considering the non-linear term involving y”, 


we shall, as a first approximation, neglect the influence 

of the harmonies and retain only those terms involving Wo ° 
Thus .cubing equation (4), reducing terms to first 

power functions of w,t, and discarding all functions of 


AWet, 5 wpt, we obtain: 


yx” ss 3/4 a” gin wt + 3/4 de cos Wot + 3/4 a Wt 


(7) 
+ 3/4 e” cos wet 


Substituting #, y and y° from equations (5), (6), and (7) 
into equation (3) and equating the coefficients of the 


sine and cosine terms respectively, we have: 


2 
-awss - LNCWp + wn + 3/4 @, ae (a + ”) a Ww Sst 


P (8) 


~cws + 2nad, * rete + 3/4 we e“e (a” + e”) = @ 


. 
Rewriting: 


2 (wi - we) - 2nw,c * 3/ 4 wo *2 a (a” * e”) _ wo” Yst 
\e 
4 (9) 


Cc (we ~ w-) * 2nW@pa + 3/4 wi be eC (a” + o*) = O 


3 
=a 





el 


Let R= — , 


5 « on 

w, 
K 

K = 3/4 ; 


Substituting these relationships into (9), we obtain: 


a(leR”) = YRe + ax(a” + ce”) = Vet 


e(1-R”) + JRa + xia” + os") = Q 


Rearrang ing: 


2 


f{l-k +* XK (a” + e) la- JRe = Vat 


(10) 
fil-R?+K (a +e’) Je + S Ra = 0 


2 
Let S = a” +e = the square of the amplitude for steady 


atate conditions. 


Substituting and operating on (/0): 


2 


A KS)a J Re + sna = Vot 


(1 = R° + xs)" - 2(1 - RB 


(1 - n® + xs)"e” + 2(1 - R® + KS)e Y Ra + J”R@a” = © 





(1 = R~ + xs)*¥5 + Yy*n®s = yZ, (11) 


This is a cubic equation in S, the real roots of whieh will 


ape © ere teee ce 68 el el bees mee & Uf ot 
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give values of the possible steady state amplitudes. The 
criteria for stability of these amplitudes will not be dis- 
eussed here since the regions of stability and instability 
of the non-linear response curve are well known and are 
covered in most vibration texts. 

Since solution of a cubie equation is unhandy, §& 
may be made the independent variable and the equation ex- 


pressed as a quadratic in ?. 

Rewriting equation (11): 
2 
st/. 


(1 + xs)” ~ 2(1 + Ks)R® + R* + YR? ¥st/ he 
s 


2 
[(1+*#Kks)-R®°} +* Yy*R® = y 


4 


f 2 
R os 2(1 + KS « Th yRe = Fst / + (i + KS)” = Q 
5 


whence 





2 
Re = wBe(Ie dA) USE ga yeh) Ksva- veh "Way = i= ans = 1Ps* 


2 “yy, * : 
R© = KS + (1- 0%) = [Yar p= UK *(1 = 3A)°-1 (12) 


Equation (12) is the expression used in the conpue 
tation of the non-linear coupling response curves. It lends 


itself readily to a computing form. 
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3. Other theoretical considerations 
(a) Evaluation of Ty 


From the preceding derivation; 
Tilo 
"Hee 
fhe net torque exerted by the motor must be equal 
to the resisting torque of the generator plus the damping 
torque, Th and Tp respectively. This net torque is a cone 
stant at dise "a", out in being transmitted through the 
linkage of the Whitworth drive aequires an oscillatory com- 
ponent. From the mechanics of the system it may be deter-~ 


mined that 


» © p* - 2p cos (wt) + 1 
bd a 


p™ = p cosa (wt) 


where p = = 


r = the distance of the driving pin of the linkage 
from the center of disc "a", equal to 1.5 inches 
for this installation. 

e = the eccentricity between discs "a" and "b", equal 
to 0.125 inches for all runs made. 


T." Tr * Ty 


T,' = total driving torque applied to dise “b". 


om * etn? 


ad ares eet ld 
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Evaluating the constants and expressing as a series 


we have: 


_* T,(1.00-.085 cos wt - .003 cos 2wt - .0001 cos Jwt 
+ eeweeeee) 
Disearding the harmonics above the first as negli- 
gible, there remains, superimposed on the constant torque, 


a simple sinusoidal torque of maximum amplitude: 


Since vibrations of the masses "B" may be trans- 
mitted through the Whitworth drive to masses "A", the latter 
will be caused to oscillate also. These oscillations result 
in an inertia torque being impressed on dise "b". Desig- 


nating this torque as rT, we may writes 
> fo om 
a hs 


Again by the mechanics of the “Yhitworth drive, sub- 


stituting the constants used above, we have: 


X , =m" (.083 sin wt + .007 sin 2wt + .0005 sin 3 wt 


we ot teeees) 


Neglecting the harmonics higher than the first and 


substituting above we obtain: 


T,' = .083 m," I, sin wt 
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Or, designating the maximum amplitude of the sinus- 
oidal torque as T,, and noting that wis the same as the 
speed of rotation of the coupling (vin., Wp), we have: 


e 


The oscillatory component of the torque applied 
to system "1" will be the sum of the oscillatory component 
of the driving torque plus the inertia torque. This will 


be a sinusoidal torque of maximum amplitude: 


Ty . Ty + Ty 


Substituting in the expression for T,;, we have: 


. + +7 te 
tT, = [083 @, 1, + .088 (1, wr 


+ I, 
- [085 w,” (x, + 1,) + .083 (7, + t>)] a2 
' a ee 'y) | rep : 


Ty is evaluated from a speed vs. load-torque curve. 


T.. ls evaluated from the hysteresis loop for the par- 


oO 


ticular set-up being considered. 
(b>) Evaluation of the damping factor, X . 


The damping factorsused in the computations for 


both the shafts and the non-linear couplings were s0 deter- 


ee ne 
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mined as to make the maximum amplituce of the theoretical 
curves at resonance the same as that observed. Although 
this is apparently a subterfuge, its use will be defended 
in the discussion to follow. 

For the linear system at resonance, b= 5 and 
may ve determined by computing backwards from a given maxi- 
mum amplitude. 

For the non-linear system 5 may be determined as 


follows: 


Differentiating equation (11) with respect to f, we obtain: 


2 2 a& of ai 
| (a = n®) + xs]? SS + es [(a- w®) + Ks) (-en + KG) 
2.2 a8 a 
+ OR 2 + 2 = 
R yo RS = 0 


2 
ae £0 -R*) #2 xs (1 - KR”) + Ks” + @ Ks ( 1 - B®) 


+ 2 Ks” + x2? ]eans[ta- a) +n -& | 


or: 
as 2RS (2-2 R” + 2 ks - } 
Oe eS 
aR (l= R°)’ + y"n® + 4 5S (1 = RR“) + 3 XS 


The maximum point on the resonance curve occurs when the 
slépe is zero. Since K and 3 are not zero, it is neces- 


sary that; 





2728 
ea-e2Rr°+exns- J” = O 


1-R° + KS- 3, = 0 


Substituting for R® from equation (12): 


c 2 - £ - 
1-KS-1+ Up Ft V ye st / = B°KS+* (1-8°/o ) 73 +K8~ I" /p 0 


whence 
yee, - Vexs + (1 - 0H) -1 = 0 
S 
yx” ~ YKStie- y+ 0%, -1 = 0 
st/. 


S44 = (KS +1) 5" + yey = 0 


2 
= 2 eS 
8 2 (KS * 12 Veg + 3)®. y¥° ) (13) 











DISCUSSION 


Sinee the comparison of the non-linear with the 
linear system was made on the basis of an “equivalent shaft" 
it is necessary first to define this term. This definition 
follows: “A shaft is considered equivalent to a given 
coupling when the slope of its torque-deflection curve is 
the same as the initial slope of the torque-deflection curve 
of the given coupling.” On this basis the initial flexi- 
bility of a ceupling and its equivalent shaft is the same 
(1.@e, k of the shaft equals k, of the coupling). Thus the 
two systems are brought together at an initial point. De- 
parture of the behavior of one system from that of the 
other as operation occurs removed from this point (viz., 
amplitude of vibrations become greater than sero) must be 
due to whatever difference exists in the spring character- 
istice and the damping of the two systems, since the moments 
of inertia remain essentially the same. 

As the purpose of this investigation was to deter- 
mine the effect of the non-linear spring characteristic, 
it would have been desirable to make the damping of the two 
systems identical. However, no attempt was made to do this, 


ang it will therefore be necessary to estimate the effect 
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of the difference in damping in order to draw a conclusion 
as to the effect due solely to the non-linearity of the 
spring. It is doubtful that a coupling could have been de- 
signed such that the damping in the coupling system and in 
the equivalent shaft system would have been identical. In 
neither case was the damping purely viscous; rather it was 
of some complex type, evaluation of which would have been 
quite difficult. Following aceepted practice, the damping 
was expressed as an equivalent viscous damping, in amount 
equal to that necessary to give the observed maximum ampli- 
tude at resonance. The method of making this evaluation 

for each system has been explained under "theoretical con- 
siderations.” It is believed that the major cause of dis- 
crepancy between theoretical and observed values of ampli- 
tude is due to this approximate means of evaluating the damp- 
ing, anc it is therefore suggested thet if future investiga~ 
tions are undertaken, they be directed primarily toward a 
more aceurate solution of this problem. 

Accepting the demping determined as described above, 
it was still necessary to arrive at some basis of comparison 
whieh would present the true effect of the non-linear spring 
characteristic alone. Since the damping in the coupling was 
much greater than that in the shaft, it is obvious that s 
comparison neglecting the difference in damping wonld present 


the former in too favorable a light. Therefore a theoretical 
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maximum amplitude was computed using the values of damping 
ana torque at resonance of the equivalent shaft. The maxi- 
mum amplitude thue obtained was almost double that observed, 
but still about half the maximum amplitude of the equivalent 
shaft, aa may be seen in Figures 6 to 9 inclusive. 

Damping also complicated the evaluation of the 
torque-deflection curve of the non-linear system. Since 
the deflection was changed very slowly in obtaining the 
curve, the only damping present was that independent of the 
speed, viz., Coulomb damping. The result of its presence, 
however, was to give a hysteresis loop rather than a single 
curve. Inepection of the curve indicates thet the friction 
was greater with decreasing torqne than with increasing 
torque. Thus we have the spring characteristic confused 
with part of the damping. The true spring characteristic 
must lie somewhere within the loop. A characteristic evalu- 
ated to pass through the maximum point of the hysteresis 
loop was found to be unsatisfactory since, at the larger 
angles of twist, conditions are such as to make the damping 
efféct overshadow the spring effect on the curve. Similar~ 
ly, the mean line between the two boundaries of the loop was 
found to give poor results because it was unduly iowered by 
the Increased damping effect of the descending branch of the 
curve. By trial it was found thet a curvé approximately 


parallel to the upper branch of the hysteresis loop, passing 
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through the origin, gave the most uniformly satisfactory 
results. Since the initial slope of the approximate curve 
represented by the cubic equation is fixed (it is equal to 
the initial slope of the upper branch of the hysteresis 
loop), it can be foreed to coincide with the actual curve 
at only one point other than the origin. This coincidence 
was made approximately at the mid-point of the actual torque- 
deflection curve. Since the actual curve itself is at best 
only an estimate as shown above, a small departure from it 
by the approximate curve is not considered serious. 
Referring to equation (12), it should be noted that 


7 dis a function of W {or of WH), which is in turn a func~ 


E 
tion of R. The equation is thus not in proper form for 
evaluation as a quadratic. ltowever, a value of i] may be 
estimated, from which T, may bé@ computed. If the value of 
R thus determined gives an N which differs by more than 10 
per cent from the N estimated, the computed may be used 
for the next estimate and the value of KR recomputed. If 
there is no other means of estimating KN for the first ap- 
proximation, the value of N, (which is kmown for any par- 
ticular system) may be used for ell but the lowest ampli- 
tudes. If the first estimate of N is considerably lower 
than the actual value, the quantity under the radical be- 
comes negative, giving a eompliex number for R°. In general, 


knowing the value of N,» and having a knowledge of the gen-~ 





- — - ae") ° 
4 ’ ‘( ™~e = ° 
"en ) 6 . 
‘ : - 2 
—s tie | 
. —— 
; 
= > i 
— ~ ae a | 7 
Ps a — 7 =. rv a oo al 
. CS 
a a — - 
— all oa ar ee 
7 atte eee “ilies asi an an a 
7 
1 as 7S —_— «4 
i — Sa ee ad Te 










J tre pany 


a 






wer 
“at oa 
AP pte mt tee 
ot = 0 Ree es et 6 pe et eee 
2 © eet + eet et yt © eer eee 









es 





35 


eral shape of non-linear response curves, a sufficiently 
satisfactory estimate can be made so that not more than two 
approximations are necessary. A 10 per cent variation in 
the estimated value of H will sive approximately a i per 
cent variation in the computed value of N. From this it is 
apparent that the result is net being forced by the esti- 
mate. If equation (//) had been solved as a eubie in S it 
would not have been necessary to make the above approximae 
tions since in that case FR would be the independent variable. 
It 1s considered, rowever, that the above method lends it- 
self more reedily to masse calculation (ieee, a computing 
form). 

The observed amplitudes indicated on the response 
curves, Figures 6 to @ inclusive, were obtained directly 
from the readings of dials "b" and “ad” and represent the 
angle of twist from the mean position, viz., half the sum 
of the maximum angles of twist in each direction. 

It is interesting to note that in the cases where 
the non-linearity is most pronounced, the overlap of the 
upper and lower branches of the response curves is suffi- 
cient to allow advantage to be taken of the effect, as fole 
lows: If it is desired to run at a speed approaching the 
eriticail, increase the speed first until the critical is 
passed, then reduce the speed to that desired. Thus the 
machine may be made to operate on the lower branch of the 


curve where the amplitude is much smaller. 
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CONCLUSIONS 


The results of this investigation indicate that a 
non-linear spring characteristic in a torsional system is 
effective in reducing the maximum amplitude of vibrations. 
An exact determination of the mapnitude of this reduction 
was complicated by the existence of complex damping in the 
systems compared. However, s reasonable approximation seems 
to be that for systema of the type investigated the non- 
linearity will cause a reduction of about 50 per cent in 
the maximum amplitude from that of the comparable lineer 
system. 

It was also demonstrated that in a non-linear sys- 
tem advantage may be taken of the overlap of the upper and 
lower branches of the response curves to further reduce the 
maximum amplitude which must be accepted when operating in 
the critical zone. For the cases considered this additional 
reduction amounted to about 2O per cent of the maximum ampli- 
tude indicated by the non-linear response curve. Keferring 
to Figures 6 to 9 inclusive, it is seen that the response 
curve of the non-linear system shows apprectable amplitude 
over a wider range of frequencies than does that of the 
linear system. However, the great reduction effected in 
peak amplitude is considered to more than compensate for 


this disadvantare in most installations. 
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APPENDIX I 


Sample Calculations 
In the following we shall consider first some gen- 


eral calculations used for all computationa. Following 
this, sample calculations for the 4-1/5 inch eoupling are 
presented, together with those for its equivalent shaft. 
Finally, computations are made for the theoretical maxi- 
mum amplitude considering the shaft torque and damping to 
be applied to the coupling system. 


General Calculations 
= J, * I. + i, + I 


ts c 


= .01709 + .02650 + .023554 + .01085 


z ,078 slug ims” 


= .01085 + .06594 * ,00125 


= .O78 slug in.” 


SE 





o7 


I, ig 


+ 
4 Tp 


I; 


2078 x .078 2 
078 + 078 e059 slug in. 


ae = os. 
2 I+f .078 + .078 2 


Calculations for the 4-1/3 inch coupling 





Evaluation of the approximate torque-deflection curve. 


From Pigure 11, the initial slope of the torque- 


GQeflection curve is: 
k, * «534 in./o x 57-3 = 50.4% in./rad. 
and T = 17.6 in. when ¢ = 20°, 
20 


> * 57.3 = 349 radians 


ic 
2 3 
. * ky ° ky, 7 ) 
3 
17.6 = 3004 (2346 + oe x as} 
i 


= 5045 


Bes 


4 | 
K = 3/4 K « 3/4 x 5.43 = 4.1 





Determination of the natural frequency. 


o,- Vk] =~ V30.4 = 28 rad./sec. 


@. * 265 PeDel%e 


Kvaluation of the damping. 
The maximum amplitude observed for this coupling 
wae 29.5° at § = 387. 


= n= Re = 40.5 rad./sec. 
= .085 [ w, (1. + I.) * (7, + J 2 
7”? Se Sy 8 tS 
: 


= .083 [20.5 (.01709 + .02650) * (2.50 * 1.20) Ix 


- 3.154 ine 


Tr oelb 


of « 3o35 2 .1036 
Yet ik, 3004 E 
Y,. 7 +0107 


2960 
O7 ed 


2.7r = 1264 


= 0014 


9 
a 


w 
max Ymax 


2 
¥st — 00107 = 4.0404 


Pay ee 


“max 
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y? = 2 (KS +14 V (KS + 1)? - vats. ) (18) 


Substituting: 
x2 = 2 (4,1 x 264+ 1+ V(4.1 x .264 + 1)" - 0404 ) 
x” = .0196 
Y= .1400 
Evaluating the response curve: 


ere Eee ee 
R® = KS + (1 - 3%) + V Yst/, ~ ORS * (1 - S/o)" - 1 (12) 


Substituting: 
R° = 4.18 + .990 + Y 7.” - .088 = .020 
2 en en 
9238 


ve may now substitute values for S and solve for eg. 


For convenience, take the amplitude y = .1 radian (é = 5.73). 
Then S$ = y° = .01 

To evaluate Ty it is necessary to estimate a value of NN. 

A substantiation of this procedure is given in the discussion. 
For y = .l1 , for the upper branch of the response curve 


N = 220 RPM 
est 


For the lower branch of the response curve 


N = 370 RPM 
est 
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Preceeding for the upper branch of the response curve: 


Ty 


1” 


T, = .083 [we (I, + I.) (tT, + T)) 1 


= .083 | (220m) (.0436) # 2.76 ] x3 

# 1.09 
T, and Ty were ebtained from Figs. 18 and 11 respectively. 
Substituting in the equation for R® and using the negative 


sign before the radical since the lower of the two values 


for R® 4s desired: 


2 [1.08 
r e * @ ——_ - 6 s => «@ 
R 4.1 x .O1 990 955 = TOT 08 x .01 020 


a 


1.031 - .3289 
=  .7021 

R= /.YWOST = .837 

x = RM = ,837 x 268 = 225 RPM 


Proceeding similarly for the lower branch of the response 
curve, and using the positive sign before the radical, we 
ovtain: 

R= 374 


Thus tio points on the response curve have been obtained, 
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Vize? 
QO 
P= 5.73, N = 225 
B= 5°73, n= 374 


Following the same procedure sufficient points may be obe 
tained to determine the theoretical response curve for the 


4-1/3 inch coupling given in Fig. 7. 


Computations for the equivalent shaft. 
Using calculations given above for the coupling: 


#4N/ a 


k ® Ky = $0.4 
W, = 27.8 rad/sec 
Cm = 268 
At resonance; 
268 77.” 
T, = 083 | (228-0) (.0436) + 5.08 | x1/e= 1.54 
P ove. = 100° 


y = 1. 1.747 rade 


D708 


ed 
i 
® 
“a 
¥-94 
— 


maxX 
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1.54 
30.49 


1.747 = 
¥ = .0292 
Values for the response curve may now be computed: 
Determine the armlitude at WN = 250: 


250 as é 
= = od 9 R - - 0933 = 03 
Ty = 1655 San Fe 70 


as SNe 


me Va ht oe 
(1 - n°)? + yp 


1655 /30.4 


(i - 237) + OO98 x 870 


2 ,554 rads. 
0 
P= .334 x 5765 = 1901 


fhis gives one point on the linear response curve of Fig. 7. 
Other points may be determined similarly. 

Computations for the theoretical maximum amplitude 
of the non-linear system when the torque and damping at 
resonance of its equivalent linear system is applied to it. 

From previous computations - 


¥ = 0292 


Under “theoretical considerations it was shown that 


ad 


od 


'* an an 2 
ak 

oe am 

~~ a nee - 
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el 
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- e 
—_~el 


‘oat’ bemwe cee °° coe? coke et! : 


“ze 


at maximum amplitude of the non-linear res»onse curve < 
4 2 2 
8/4, - (KS +1) 9 * yst/, = 0 

Substituting and solving for 5 - 


— 
.opon* - (4.18 ¢ 1) x .0292 + 1.54 = 0 


a> 


S = 728 
y = YS == .855 radians 


f= .853 x 57.3 = 48.9 See Figs 7 
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Data 
Koments of inertia; 
I, ~ 01709 slug in”, 
Ig = solgs * * 
I =z  .02650 - . 
I = .02354 ss . 
._- —_ * * 


I. = ,06594 "* 7 


"Kquivalent" shafts: 


All “equivalent” shafts were of mild steel, 3/16 


inches diameter. The lengths were as follows: 


58 inch shaft equivalent to 4 inch coupling 
45-1/2 Ht tt " +# 4-1/3 tt tt 
50-5/8 +=" - : *“-2/s * . 


60-7/8 t "n tt " 5 f 1% 





TABULATIGCN GF CONDERGED DITA ARD COPPRESOIDING COMPUTED VOLUES 
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“h” §ndieetes operation on lower brunch off nom-lineer response curve. 


400 








4e9 









"A* ¢ndiestes opercetion om upper branch of 


S404 1705 1603 < 14.0 14.5 
5504 19.5 17-& i oan a —s 
400!. 2220 2305 $204 Ss Oe 
4284 2520 2600 | 40a 25.5 2367 
430 2820 2603 BUOA 2668 2605 
$80 10.0 90 370". m0 28-0 
4000 70 7.0 380" 26 2901 
4500 405 60 0 | —sBA 29.5 2904 
3400 1205 9.6 
350D 8.8 8.0 
400D 6.0 5e0 
138 ineh She 1/2 Shalt 
ee Se arse roan 
imal olitude | amplitude | TEUMMEY | seoiitude | Amp’ Litude 
" -seheareas | Toiaghees [IE TG jadoprecs | got 
260 Bed Te? 200 3.8 Se 
285 13.0 12.3 240 12.0 1240 
5 27-5 18.5 250 205 19 
= 5608 37.05 250 475 55<0 
305 1665 3205 285 18.0 1903 
325 11.0 15.6 | 500 12.5 13.3 
370 2.6 %.5 . 400 $.0 500) 
400 560 Be? 


norlinear response Curves 
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